program hydrostatics

    implicit none
    real :: rho, g

    rho = 1000.0
    g = 9.806
    print *, "p = \rho*g*h =", rho*g*1.0, "Pa"

end program hydrostatics
!> p = \rho*g*h =   9806.00000     Pa